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ABSTRACT 

We investigate the sources and amount of dust in early galaxies. We discuss dust 
nucleation in stellar atmospheres using published extended atmosphere models, stel- 
lar evolution tracks and nucleation conditions. The Thermally Pulsating Asymptotic 
Giant Branch (TPAGB) phase of intermediate mass stars is likely to be the most 
promising site for dust formation in stellar winds. We present an elementary model 
including dust formation timescales in which the amount of dust in the Interstellar 
Medium is governed by chemical evolution. The implications of the model for high 
redshift galaxies are investigated and we show there is no difhculty in producing dusty 
galaxies at redshifts above 5 if supernovae are a dominant source of interstellar dust. 
If dust does not condense efficiently in SNe then significant dust masses can only be 
generated at z > 5 by galaxies with a high star formation efficiency. This is consistent 
with the high SFR's implied by sub-millimetre sources found in deep SCUBA surveys. 
We find the visual optical depth for individual star forming clouds can reach values 
greater than 1 at very low metallicity (1/100 solar) provided that the mass-radius ex- 
ponent of molecular clouds is less than two. Most of the radiation from star formation 
will emerge at IR wavelengths in the early universe provided that dust is present. The 
(patchy) visual optical depth through a typical early galaxy will however, remain less 
than 1 on average until a metallicity of 1/10 solar is reached. 

Key words: ISM:abundanceS"dust,extinction-galaxies:high-redshift 
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1 INTRODUCTION 

The role of dust in galaxies at present and at earlier epochs 
is still highly controversial. Recent dust models have at- 
tempted to understand the evolutionary cycle of dust (Ed- 
munds 2001 (hereafter EOl); Dwek 1998 (hereafter Dwek98); 
Tielens 1998) but there can be no real advancement whilst 
there still exists inherent uncertainties over dust sources, 
composition, destruction and the dependence of dust pro- 
duction on metallicity evolution. It is well known that dust 
production in the Interstellar Medium (ISM) is too slow 
(timescales for dust growth are smaller than destructive 
timescales (Tielens 1998)) so formation sites during stellar 
evolution are required to explain the vast amounts of dust 
we observe in the universe. 

It is generally accepted that dust is formed in the atmo- 
spheres of evolved stars and ejected via stellar winds with 
mass loss rates of up to MQyr~^ (Gehrz 1989; Whittet 
1992) and contribute up to O.SM© of dust to the ISM in 
our Galaxy per year (Bode 1988). Others however, propose 
that supernovae explosions must be an important contribu- 
tor to dust formation (Dwek & Scalo 1980, Todini & Ferrara 
2001 (hereafter TFOl)) to account for the dust observed at 



high redshifts (as seen in Damped Lyman Alpha systems 
and in SCUBA data (Frayer & Scoville 1999)). The lack of 
significant XMM detection of SCUBA sources in a deep field 
(Waskett et al, 2003) imply these sources are probably not 
AGN and even if they were, then they are likely to be heavily 
obscured and at high redshifts {z > 2). Almaini et al (1999) 
also show no significant Chandra detections of other sub-mm 
fields. If not AGN then they are probably extremely dusty 
starburst galaxies, implying that dust is already present in 
the high redshift universe. It is difficult for the dust to have 
originated from low-mass evolved stars at redshifts > 5 as 
their evolutionary timescales (10* - > lO^yrs) are compara- 
ble to the age of the Universe at that time. So it is argued 
then that if dust is present in the early universe, then Type 
II SNe are responsible. However, some estimates of the con- 
tribution of SNe dust mass to the current ISM show an injec- 
tion of only 10"" A/© per year (Bode 1988; Gehrz 1989), and 
a recent study of the InfraRed emission from three supernova 
remnants (SNR) has shown no evidence of dust formation in 
SN ejecta - in fact, it was found the emission was probably 
circumstellar in origin (Douvion et al 2001a). 

Even dust production in stars has its own complications. 
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Depending on the star's composition (whether it is carbon 
rich (C stars) or oxygen rich (M stars) and what fraction 
of a stellar population is in these stars, there will be a dif- 
ference in the returned fraction of carbon and silicate type 
grains to the ISM. The problem compounds when O rich 
stars may become C rich as they ascend the Asymptotic Gi- 
ant Branch (AGB) whilst undergoing pulsations (Marigo et 
al 1996). IRAS observations (Molster et al 2001) and obser- 
vations of AGB stars in the Large Magellanic Cloud (Trams 
et al 1999) show a high degree of silicate dust in C stars and 
it is yet unclear whether this is due to silicate formation in a 
previously O rich star or due to an O rich binary companion 
(Molster et al 2001). 

The questions posed and discussed here are then: when are 
the conditions for dust formation reached in stars? What is 
the source of dust at high redshift and is dust production 
dependent on motallicity? How much dust is present in these 
early galaxies and how will this affect observations? The aim 
of this paper is to ask whether or not we can build up a sig- 
nificant amount of dust in high redshift galaxies using the 
available literature and models. 

In Section 2 we use published stellar atmosphere and 
evolution models to determine when dust may form in stars. 
We then estimate the efficiency of dust formation in stellar 
winds and supernovae. Section 3 presents the well estab- 
lished chemical evolution equations and their applications 
for following dust evolution in a galaxy. In Section 4 wo at- 
tempt to understand the source of dust in early galaxies, 
using the modified elementary model from EOl. We incor- 
porate the finite lifetime of stars into the dust model and 
compare this with evolutionary timescales to determine the 
maximum redshift after which dust can be observed in the 
Universe for different star formation rates and cosmological 
parameters. Finally, Section 6 gives the conclusions obtained 
in this report and an outline of future work. 



2 DUST PRODUCTION IN STARS AND SNE 

The majority of work has assumed that the formation of 
solids from a gas in stars can be described using classical 
nucleation theory. Condensation of a species X occurs when 
its partial pressure in the gas exceeds its vapour pressure 
in condensed phase (Whittet 1992), with subsequent parti- 
cle growth by random encounters leading to cluster forma- 
tion. Chemical reactions will further assist cluster produc- 
tion (Salpeter 1974). The rate of grain growth will there- 
fore depend on the temperature, Tg and pressure, Pg of 
the gas and the condensation temperature of a species, Tc- 
The optimum conditions for nucleation are thought to oc- 
cur in the range 10"* < Pg{Pa) < 10^ and Tg < imQK 
(Salpeter 1974; Yorko 1988). This gives very strict restric- 
tions on where nucleation can occur and this is thought to be 
satisfied in giant atmospheres at certain stages of evolution. 

The type of dust formed will depend on the relative 
abundance of carbon or oxygen with carbon dust able to 
form in temperatures Tg < 1800-/^. For silicates to form via 
condensation we require temperatures below 1500K (Bode 
1987; Whittet 1992; Salpeter 1974). These conditions repre- 
sent the boundary 'nucleation' in Figures 1(a) and (b). 



2.1 Nucleation in Stellar Atmospheres 

Unfortunately there is not yet nearly enough hard observa- 
tional evidence to follow through the whole process of dust 
formation in stars so, in an attempt to investigate when 
nucleation conditions are reached in stars, we examine ex- 
tended atmosphere models for K and M giants (Plez et al 
1992a & 1992b; Hauschildt et al 1999b; Brown et al 1993). 
We use these models to estimate the temperature and pres- 
sure in the gas out to a stellar radius where the contin- 
uum optical depth at 500nm has fallen to 10~®. The at- 
mospheric conditions depend on the mass of the star, M, 
its surface gravity {g = ^^), and its effective temperature 
at the surface, Teg. The published models assume spheri- 
cal geometry with homogeneous stationary layers in local 
thermodynamic equilibrium (LTE) and use numerical codes 
to simulate the gas properties in stellar atmospheres (e.g. 
PHOENIX (Hauschildt et al 1999a) with these models avail- 
able from the website: http://phoenix.physast.uga.edu/NG- 
giant). Although this is probably not a realistic description 
of conditions in stellar atmospheres where pulsations, mass 
loss and inhomogeneities will surely exist, they will serve to 
make an estimate of the mean condition for the processes 
occuring. 

Model atmospheres for C stars were obtained from 
Jorgensen et al (1992) and an observational model (Harper 
et al 2001) of the red giant AlphaOri is included for com- 
parison (with Teff = 3600K, M = 12Mq, log{g) = 0.0). 
The model P-T curves are shown in Figure 1(b) for Teff > 
3000K, and Figure 1(a) for Te^ < 30Q0K. They clearly show 
the models do not extend to the low gas temperatures ob- 
served in AlphaOri, and we will assume the boundaries do in 
fact extend downwards in a similar manner. It is also clear 
that as Teff increases, the likelihood of the gas temperatures 
reaching the nucleation boundary decreases. The figures in 
Hauschildt et al (1999b) show for various values of log{g), 
that only those atmospheres with T^g < 4000fc will reach 
gas temperatures below 2000K and hence attain C nucle- 
ation conditions. In order to reach the silicate nucleation 
boundary, we require T^ff- < SOOOA' (Fig 1(c)). 

To determine to which stars these model atmospheres 
apply, the luminosity was calculated using Equation 1 and 
simply placed on a Herstzprung- Russell diagram using their 
respective T,,ff values (Fig 2). 

L. = 4naT^ff (1) 

The stars do indeed correspond to the evolved stages 
expected, but for dust nucleation properties to be reached 

we require the stars to lie on or just beyond a late AGB 
phase (evolutionary tracks of solar metallicity, Zq = 0.02, 
from Schaller et al (1992)). It appears then that the stars 
seen to produce dust using the arguments given here, only 
satisfy the conditions at the very latest stages of their 
evolution which is only a very small fraction of their life 
(lO'' - lO'^years). This phase is the TPAGB (Thermally 
Pulsating Asymptotic Giant Branch) as described further 
in Section 2.2. The stellar evolution models from the the 
Padova group (http:/ /pleiadi. pd.astro.it) also give the same 
conclusion. The stellar tracks can reach lower cfi^cctivc tem- 
peratures (Pols et al 1995) depending on the opacitites, 
overshooting and mixing parameters used although a recent 
study of these models suggest there is still more to be done 
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Figure 1. The gas properties of model atmospheres of different stars from available literature models (see text). The RG/AGB curves 

represent the temperatures and pressures seen in the atmospheres of evolved giants if (a) T^g < 3000K and (b) T^g = 3600 — SSOOX. 
The solid curve labelled solar represents the atmospheric properties for a IMq star with logT^g = S.SX (Hauschildt ct al 1999b). An 
empirical atmosphere (Harper ct al 2001) representing the giant star Alpha Ori and the nuclcation boundary of carbon are included in 
both figures. For a higher effective temperature, the ability to reach the nucleation conditions in stellar atmospheres decreases (b). The 
nucleation conditions in a carbon-rich gas as adapted from Yorke (1988) and Salpeter (1974) are shown in (c). For an O-rich environment, 
silicates can be formed if Tg < ISOO/sT for the same pressure range. The arrows represent the direction in which nucleation conditions 
are satisfied for the elements shown. 



in this area (Young ct al 2001). Wo await further progress 
on new tracks given by this group using their new stellar 
evolution code, TYCHO. 

In further sections wc assume that only those stars with 
Teff < 4000i^ and which lie on this evolutionary branch 
are significant producers of dust (although this is actually 
a maximum limit since these are the conditions for carbon 
nucleation and we would expect to require much lower tem- 
peratures for the formation of silicate dust). 

It is important to note that because of this reasoning 
we have ignored dust production in the atmospheres of stars 
with Mi > 25Mq. ObservationaJly however, dust shells are 
seen around a class of massive stars which suffer large mass 



loss lO~^M0jyr~^) in the late stages of their evolution. 
These arc the Wolf-Rayet stars, which have three phases 
characterised by the dominant species of O, C and N in 
their spectra. Dust shells are only seen around the very cool, 
carbon rich WR (Dwek98), and will therefore contribute to 
the production of carbon dust in the ISM. Could this be 
a significant dust formation site? Recent observations of a 
WR binary system have shown that the formation of dust 
in these stars may be increased via wind-wind interactions 
with at least 20 per cent escaping and thus enriching the 
ISM (Marchenko ct al 2002). Their rcsuhs find a total dust 
mass of 2.8 x 1O~^M0 which coupled with the mass loss rate 
gives an estimate of the condensation efficiency of dust in 
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Figure 2. The HR diagram with added stars corresponding to the model atmospheres used. The triangles represent the K & M giants 
whilst the crosses represent the C stars. 



WR of ^ 6 per cent. But WR stars are also very rare in 
galaxies as compared to low-intermediate mass stars, with 

the fraction of stars in a given population with M > 25Mq 
^ 10^'' implying that the yield of dust from such a process 
will not be significant. 

2.2 Mass Loss and the TPAGB Phase 

The dominant processes in stellar mass loss are pulsa- 
tional and line driven winds with mass loss rates ^ of 
- 10-"* - IQ-'^MQyr-^ and 10"'' - IQ-^^M^yr'^. The im- 
portant question to ask is how much of these outfiows are 
in the form of condensed material? Whether or not dust is 
formed from the gas will depend on the temperature and the 
density inside the wind (Section 2.1). A rough estimate of 
the dust production in these winds can be made by investi- 
gating condensation conditions in stellar atmosphere models 
which in turn are tied to stellar evolutionary models which 
predict Teff and M throughout the stars evolution. 

The mass lost throughout the star's lifetime must sat- 
isfy the initial-final mass relationship shown in Weidermann 
(1987) and Marigo (2000). Thus for low-intermediate mass 
stars 1 < Mi < 7 Mr.), the mass lost in stellar winds and 
novae is the difference between the initial mass, Mi and the 
final remnant mass. Ma. For high mass stars (7 < Mi < 
40Mq ) the mass is lost via both winds and SNe. Stars with 
mass beyond 40Mq may become a black hole and swallow 
up the elements produced (Maeder 1992) so they are ignored 
here. 

We have followed the evolutionary model grids in 



t Dust driven winds only become important at high metallicities 
where there is sufficient dust production to drive the gas outwards 
(Willson 2000). 



Schaller et al (1992)'f, up to the end of central C burning for 
high mass stars, the early AGB phase for intermediate stars 
and up to the He fiash for low mass stars. Some of the mass 
lost during a stars lifetime must occur during a very late 
AGB phase or after (the mass lost from the published mod- 
els does not explain the difference between the initial mass 
and remnant mass§. In the TPAGB & stellar wind phases, 
M &i T were taken from Blocker (1995a) and Marigo (1996), 
but how and where does this occur? As a star ascends the 
AGB its thought that (depending on the previous mass loss 
history), it may undergo some Ho flashes and will pulsate pe- 
riodically. This is the TPAGB phase and will create a dredge 
up of material within the star to its surface. There will be not 
only a change in the composition in the atmosphere but also 
the composition of the dust (Garcia-Lario &i Pera Calderon 
2001; Herwig ot al 1999). The TPAGB phase has mass loss 
rates comparable to the AGB phase and lasts for approxi- 
mately 10'* — lO'^years depending on the initial stellar mass 
(7-3M0 respectively) and the number of thermal pulses in 
the phase (Blocker 1995a). So there is an opportunity here 
to create and inject dust into the ISM, but arc the conditions 
right? The conditions for dust formation do appear to be sat- 
isfied in the TPAGB phase as shown in Figure 3 where the 
evolutionary tracks of stars undergoing thermal pulses from 
Marigo (1996) are added to those in Schaller et al (1992) 
for low-intermediate mass stars. It may bo seen that the 
TPAGB tracks pass through the area where dust forming 
atmopsheres exist on Figure 2. The pulses may allow for pe- 
riodic dust formation by increasing the density around the 
condensation point (Blocker 1995b) and carbon stars evolv- 



t as distributed by the astronomical data center at NASA God- 
dard Space Flight Center 

§ the points of greatest mass loss up until the end of the evolu- 
tionary grids was multiplied by the duration, r, to give the total 
mass lost 
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Figure 3. The evolution of low-intermediate mass stars up to the end of C burning (solid line) and up through TPAGB phase (dashed 
line) for Z = Zq. The PAGB tracks are given in Blocker (1995b), and are not included here. 



ing on these tracks can reach T^ff < 2500Ar (Marigo 1996), 
increasing the likelihood of nucleation. The future evolu- 
tion of the star may depend on where the thermal pulses 
occur as they ascend the AGB. If the TP occurs just be- 
fore the AGB ends, the star will experience dredge up and 
may become C rich. A superwind is then thought to oc- 
cur with M ~ lO'^Moyr"^ (Schonberner & Steffan 2001; 
Blocker 1995b). If the TP occurs after the AGB ends (char- 
acterised by an increase in T^ff and a decrease in M (Fujii et 
al 2001)), then it will evolve on a constant luminosity track 
for ~ lO^years to a Planetary Nebula stage (PN) - and the 
star is now PAGB (Post AGB). This is characterised by an 
increasing Tejf and a decreased mass loss with temperatures 
reaching above 30,000K. The central star begins to ionize 
the surrounding material to create a PPN (Proto-Planetary 
Nebula) . This phase does not seem suitable for dust format 
tion. It does not satisfy the low temperatures or the high 
mass loss rates needed, so we can ignore this as a significant 
phase of dust production (although circumstellar dust shells 
are seen around PAGB stars, they are thought to have sur- 
vived from their previous AGB phase (Lorenz-Martins et al 
2001). If the TP occurs when He burning has ceased then 
the star will evolve on a white dwarf cooling track and will 
again not be suitable for dust conditions. 

So in these later evolutionary stages we have a clear 
candidate for a dust formation site - stars in the TPAGB 
and the superwind phases. Those stars which have had their 
composition changed from O rich to C rich are now more 
likely to form C rich dust. It is interesting to note that from 
Section 2.1, not only is it difficult to reach the low temper- 
atures in stars for silicates to form but we are now faced 
with the problem that those O rich stars are likely to be- 
come C rich in the most significant dust formation stages! 
This would infer that silicate dust may not be the domi- 
nant component of returned dust from stars as previously 
thought (Bode 1988), but instead SNe may be responsible 
for most of the silicate dust we see in the ISM. Dust for- 



mation in the TPAGB phase will also be highly dependant 
on metallicity as dredge ups are more efficient at low Z, 
hence the TPAGB is reached earlier on in the stars' life- 
time and more C stars are formed (Dwek98; Groenewegen 
1999; Marigo 1998; Marigo et al 1999). The enhanced con- 
version of O stars to C stars at low Z may then imply a large 
abundance of carbon dust in the early universe (although at 
Z=0.001 we need at least 1 x lO^yrs for M < 5Mq to reach 
the TPAGB phase). 

It is worth mentioning that if it is mainly carbon dust 
forming in stars during the TPAGB phase then it may be 
formed from freshly synthesised elements in the star and 
may be considered differently to dust formed from pre- 
existing heavy elements in stellar atmospheres. This may 
have consequences in models which separate the formation 
of dust into two different modes. Prom EOl, we can define 
a dust formation efficiency (i.e. how much of the available 
heavy elements will be incorporated into dust) depending 
on whether those heavy elements are from freshly synthe- 
sised (e.g. in SNe) or pre-exisitng material (e.g. in stellar 
winds) . We label these dust condensation efficiencies xi a-nd 
X2 respectively. The TPAGB phase may thus be considered 
as a xi source instead of a X2 as although the dust is be- 
ing formed in a stellar wind, its source of heavy elements is 
dredged up from freshly synthesised carbon from the core. 
Although the ejecta is likely to be mixed, the carbon dredged 
up during the TP- AGB phase is a way of producing fresh ele- 
ments at low metallicity, i.e. the production of both silicate 
and carbon dust may depend on metallicity. The implica- 
tions of this is investigated in Section 3. 

2.3 Estimating the Condensation Efficiency in 

Winds 

To estimate the condensation efficiency in winds we can use 
the models in Schaller et al (1992) and Marigo et al (1996) 
to give the mass fraction lost in winds (as listed in Table 1), 
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AMsw 


mpc 


mpz 






Total M^ust 


Implied X2 


1 




0.5 


0.5 


2.4E-5 


2.4E-5 


8.0E-6 


5.2E-6 


1.3E-5 


0.55 


3 


350 


2.4 


2.4 


3.5E-2 


3.0E-2 


6.0E-4 


1.9E-2 


2.0E-2 


0.65 


4 


160 


3.3 


3.3 


1.5E-2 


2.0E-2 


2.3E-4 


1.5E-2 


1.5E-2 


0.76 


5 


94 


4.2 


4.2 


5.0E-3 


1.215-2 


1.7E-4 


9.0E-3 


9.2E-3 


0.76 


9 


26 


7.8 


2.8 


2.7E-2 


1.7E-1 


l.OE-4 


<6.0E-4 


<7.0E-4 


< 0.7 


12 


16 


10.6 


2.6 


7.1E-2 


6.9E-1 


1.4E-4 


<6.6E-4 


<8.0E-4 


< 0.8 


15 


11 


13.3 


3.3 


1.4E-1 


1.32 


2.6E-4 


<3.4E-4 


<4.8E-4 


< 0.6 


20 


8 


17.8 


5.0 


2.2E-1 


2.73 


5.6E-4 


<2.4E-4 


<8.0E-4 


< 0.8 


25 


7 


22.3 


9.4 


3.0E-1 


4.48 


/ 


/ 


/ 


/ 


40 


4 


38.2 


30.0 


4.88 


8.01 


/ 


/ 


/ 


/ 



Table 1. The maximum mass Af^j^^t lost in winds from stars, Afg, which satisfy the conditions for dust nuclcation, for Z = Zq. The 
columns arc (1) Initial mass. (2) Main sequence lifetime in Myrs (Schaerer et al, 1993). (3) Total cjecta mass lost including SW's and 
SNc (Maeder, 1992). (4) Total mass lost in SW's. (5) Total carbon yield for low mass stars taken from Marigo (2000) and for high mass 
stars from Maeder (1992). (6) As (5) for heavy elements, Z. (7) Mass of dust from winds in the RG/E-AGB phase. Mass loss rates and 
effective temperatures used to estimate the maximum dust mass that could form were taken from the models in Schallcr et al (1992). (8) 
Maximum dust mass in the TP- AGB/superwind phase. Carbon dust formed in TPAGB assumes ~ 80% of the freshly synthesised carbon 
satisfies the nucleation conditions. (9) Total dust mass produced. (10) Implied condensation efficiency calculated from Equation 2. 



and the fraction which satisfies the conditions for nucleation 
given in Section 2.1. Equation 2 will then give an estimate 
of X2 (although this is actually a maximum as we assume 
100% condensation wherever the nucleation conditions are 
satisfied) . 

Mass into dust in stars 

X2 = (2) 

Mass of heavy elements lost in Stellar Winds 

It must be emphasised that a more detailed definition of 
the dust condensation efficiency in stellar winds would be to 
account for the free carbon in the cjecta (i.e. C not locked 
up in CO) and any available Si rather than just take the 
mass of heavy elements ejected as was done here. This has 
been done in a much more detailed manner in Dwek98 but in 
high redshift gala^Kies, it would be difficult to determine the 
composition from the sub-mm fiuxes only. Thus we try and 
estimate the maximum dust mass which could be created 
in winds and ignore detailed composition calculations. The 
new dust formed in the TP-AGB phase is likely to be carbon 
based whereas in the superwind phase, the dust is likely to 
be silicate based depending on the C/O ratio in the star. 
The majority of any pre-exisiting dust in the ejecta is likely 
to be silicates. 

To determine the average condensation efficiency over a 
stellar population we need to sum over the stellar mass spec- 
trum, with Salpeter initial mass function, </>(M) oc M~^'^^. 
The uncertainty of whether or not the TPAGB phase repre- 
sents X2 or xi does not allow a clear estimate to be made, 
and initially wc estimate X2 for the models in Schaller et 
al (1992), and X2 which also includes the TPAGB and the 
superwind phase. 

So for solar metallicity and using Table 1, 

X2 ~ 0.16 
X2 ~ 0.45 

If the condensation efficiency of nucleation < 100% then this 
would decrease x- Published estimates are 0.3 — 0.7 (Do- 
minik et al, 1993) for condensation of Si, Fc, Mg, O and ~ 0.3 
for carbon (Cadwell et al, 1994) in stellar atmospheres. The 
numbers given in Table 1 are a maximum and could be as 
low as 0.1 if we assume only 30 per cent nucleation efficiency 



when conditions arc right. This is a minimum limit for an el- 
ementary model (EOl) to reach the observed dust masses in 
the solar neighbourhood. Prom these simple arguements, it 
appears that the TPAGB and superwind phases may be the 
most significant dust production phase in the stars lifetime. 



2.4 Estimating the Condensation Efficiency in 
SNe 

Recent work on the nucleation of dust grains in supernova 
gas has begun to shed some light in this area (Kozasa et al, 
1997; Clayton et al, 2001; TFOl). Using SN1987A observa- 
tions as a basis, they all predict rapid growth of dust grains 
around 500-600 days after the initial explosion (as shown 
in SN1987A by the increase in the luminosity emitted in 
the Infrared at these times). The observational evidence for 
dust condensation in SN1987A is promising (for a detailed 
discussion see Wooden, 1997) but it nevertheless has a fear 
tureless emission spectra with no Si, C or O attributes. This 
is thought to be accounted for by iron type grains which is 
also suggested by the similarities in the iron and dust fill- 
ing factors of the ejecta (Wooden, 1997). The total mass of 
dust observed in SN1987A from IR emission is ~ 
although the elemental depletion seen in the ejecta suggests 
this should be a lot higher (Dwek, 1992). This discreprancy 
can be explained if a colder population of dust grains ex- 
ists and survives in ejecta clumps (Lucy et al, 1991). These 
grains would radiate at longer wavelengths and could pro- 
duce a further of dust not detectable in IR. Future Sub- 
mm observations of SNR's could detect this cooler dust and 
if discovered may suggest condensation efficiencies of ^ 100 
per cent. There is other observational evidence to support 
this view with dust emission seen in the young SNR Cas- 
siopiea A. IR observations of Cass A indicate that IQ^'^Mq 
of dust is present (Arendt et al, 1999; Douvion et al, 2001a,b) 
which is much less than expected from theory. Two 'recently' 
observed supernovae SN1993J and SN1994I show no sign of 
dust condensation in their ejecta and unlike SN1987A, which 
began its SN phase as a blue supergiant, these progenitor 
stars were red supergiants (Wooden, 1997). It remains to be 
seen if SN1987A and CassA are special cases. 
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Theory is still a long way from reliably estimating the 
condensation efficiency of dust in SNe of different progenitor 
mass and evolutionary phases. The most acceptable theories 
begin with nucleation via chemical or kinetic rates, but still 
can not exactly match the observations of SN1987A (with 
that of Kozasa et al (1989, 1997) predicting lower gas tem- 
peratures than actually seen). Interestingly the models in 
Clayton et al (2001) allow for a greater production of car- 
bon dust even in ejecta with C/O < 1 compared to the 
more classical nucleation models (Kozasa ct al, 1989, 1997). 
Again there may be the possibility of more carbon dust pro- 
duced than silicates. We take the elementary model pub- 
lished in TFOl to estimate the condensation efficiency in 
SNe (although this is based on the nucleation theory given 
in Kozasa et al, which has boon shown to be incorrect for 
SN1987A). The model gives the dust mass produced per su- 
pernova as a function of the mass of the progenitor star and 
the metallicity. TFOl have modelled two cases of dust pro- 
duction - (a) and (b) - corresponding to low and high energy 
SNe explosions respectively. 

To estimate xi for these two cases, we require infor- 
mation on the heavy element yield from the ejecta (as we 
are interested in finding the condensation efficiency of dust 
forming from freshly syrithcsiscd elements. Equation 3). 

Mass of dust formed in ejecta , , 

Mass of freshly formed heavy elements in ejecta 

The yields were obtained from the tables given in Woosley 
& Weaver (1995), for a grid of stellar masses and metallici- 
ties Z/Zq = 0, 10"", 10"^, 1. Weighting the values over the 
stellar mass spectrum for solar metallicity, gives: 

XI (a) -0.22 

Xi(&)-^0.23 

Other theoretical models suggest efficiencies of ~ 0.1 — 
0.3 (Kozasa et al, 1991; Clayton et al, 2001) whereas obser- 
vations suggest < 0.1 (Douvion et al, 2001b). The available 
data allows for a comparison of xi values at different metal- 
licities as discussed in the next Section, although we must 
stress that there are no observations which clearly confirm 
the extent of any condensation of dust in SNe (other than 
SN1987A and Cass A), this is purely based on theoretical 
models. 



2.5 The Dependence of Dust Sources on 
Metallicity 

Now we can ask the question does the condensation effi- 
ciency as estimated here vary with the metallicity evolution 
of a galaxy? It seems obvious that there will be less dust 

formation at low Z due to higher temperatures in the stellar 
atmopsheres but can we actually quantify this relationship? 

It is useful at this stage to re-define the condensation 
efficiencies used in previous Sections: 

• Xi " efficiency of dust condensation from new heavy ele- 
ments for the (a) high energy (b) low energy SNe explosions. 

• X2 - efficiency of dust condensation from pre-existing 
heavy elements in stellar winds of stars during their 
RG/AGB phases. 

• Xi - condensation efficiency of dust from stellar winds 



of evolved stars during their TPAGB phase if we assume the 
available heavy elements are freshly synthesised in the star. 

• X2 - condensation efficiency of dust from pre-existing 
heavy elements in the stellar winds of evolved stars during 
their RG/AGB phase and their superwind/TPAGB phases 
(i-e. X2+xt)- 

Using the same arguments as given in Section 2.3, the 
condensation efficiencies in stellar winds, X2 and X2 at dif- 
ferent metallicities have been estimated. The criteria for nu- 
cleation is the same with the mass loss rates and effective 
temperatures taken from Schaller et al (1992), Schaerer et 
al (1993), and Schaller et al (1993). The published mod- 
els allowed estimates to be made for dust nucleation for 
Z = 0.001,0.008 and 0.04^. Their corresponding X2 values 
are listed below in Table 2. Xi represents the condensation 
efficiency in the TPAGB phase if it is considered as a xi 
source i.e. freshly synthesised dust production. X2 groups 
the TPAGB into the stellar wind mode. 

For xi, the same process of estimating the condensa- 
tion efficiency seen in Section 2.4 was applied to the yields 
at other metallicities from Woosley & Weaver (1995). The 
values are given in Table 2 where (a) and (b) still represent 
the low and high energy cases. 

For zero metallicity, it has been shown that the TPAGB 
phase does occur but only at Ti-g > 4500A' (Marigo et 
al 1999). Therefore the TPAGB phase at low metallicities 
is not considered here as a source of dust production. For 
Z = 0.008, the star does appear to reach the low tempera- 
tures but for a shorter period of time as compared to those 
at solar metallicity (as shown in Marigo et al (1996)). The 
problems in trying to calculate the condensation efficiency 
at different Z become apparent when trying to determine 
the mass loss rates and the time spent in the nucleation 
conditions. With lack of any concrete information it is as- 
sumed that the TPAGB properties at Z = 0.001, 0.008, 0.04 
are similar to those aX Z = 0.02 but at lower metallicities 
less dust is formed per star due to less time being spent in 
this phase. Note that this will be a maximum estimate of x- 
There is no significant correlation between the metallicity 
and the condensation efficiency in SNe so we assume xi is 
roughly constant, ~ 0.23. The approximate variation of X2, 
X2 and Xi with metallicity is given below, 

X2(^) = -162Z^ + 9Z 

xi{Z) = -676^^ -h33Z 

xl{Z) = -5807Z^ - 612Z^ + 21Z 

3 MODELLING CHEMICAL EVOLUTION 

Using the notation given in Pagcl (1997) and Henry et al 
(2000), we can follow the buildup of heavy elements over 
time in a closed galactic system. Initially we have a mass of 
gas, g, which is converted into stars, s, via a star formation 
rate - SFR, il>{t) and an IMF, ^(m). The total mass of the 
closed system will then be 

^ as distributed by the astronomical data center at NASA God- 
dard Space Flight Center 
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Table 2. The maximum stellar weighted condensation efficien- 
cies, X2, X21 X* aiid XI in stars at different metallicities, Z/Zq. 



Z/Zq 


0.0 


10-4 


10-2 


0.05 


0.2 


1.0 


2.0 


XI 





/ 


/ 


/ 





0.16 


0.10 


X2 





/ 


/ 


0.10 


0.10 


0.45 


0.21 


x\ 





/ 


/ 


0.10 


0.11 


0.22 


0.22 


XI (a) 


0.46 


0.15 


0.18 


/ 


/ 


0.22 


/ 


XI (6) 


0.27 


0.18 


0.21 


/ 


/ 


0.23 


/ 



Table 3. The integrated yields, Px of heavy elements (pz) and 
carbon (pc) from massive and intermediate mass stars at different 
metallicities, Z. 



Mass Ranges 


Z 


Vz 




Pc 




1 - hMQ°- 


0.004 


3 X 10-3 


2 


X 10" 


-3 




0.02 


5 X 10-4 


4 


X 10" 


-4 


9- I2OM0'' 


0.001 


0.018 


2 


X 10" 


-3 




0.02 


0.015 


8 


X 10- 


-3 



" Prom Marigo, 2000 
Prom Meader, 1992 



M = g + s 



Using an infall model (such as the one described in Dwek98) , 
the extra dilution of the ISM limits the metallicity reached 
and hence the dust mass will be slightly less. Edmunds & 
Ealos (1998) show that the simple model produces the max- 
imum dust yields for a given gas fraction hence we assume 
no inflows or outflows (the effects of which are discussed 
elsewhere (e.g. Koppen and Edmunds, 1999; EOl)). The gas 
mass, g will be depleted by star formation and enriched by 
the mass ejected from stars, e(t). 

dg 



dt 



= -i^{t) + e{t) 



(4) 



The evolution of heavy elements in the ISM (Zg) is governed 
by the ejection of metals by stars (e^) and those going into 
forming stars. 

d{Zg) 



dt 



(5) 



Therefore the metallicity evolution is described by Equa- 
tion 6 

dZ _ ezjt) - Z{t)e{t) 
dt q 



(6) 



Assuming that mass loss occurs suddenly at the end of stellar 
evolution, the ejected mass and the ejected heavy elements 
are (Tinsely 1980) 



e{t) 



[m — mii(m)] ip{t — Tm)(l>{m)din 



[[m - mj{(m)] Z{t - Tm) + mpz) 
X '0(t — Tm)4>{fn)dm 



(7) 



(8) 



where mpz is the stellar yield of heavy elements from a 
star of initial mass m (values taken from Maeder (1992) 



and Marigo (2000) for massive stars and IMS respectively). 
mj{(m) is the remnant mass, Tm is the lifetime of a star of 
mass m such that a star formed at t — Tm has died at time t, 
and rriTm is the stellar mass whose Tm corresponds to the age 
of the system. Thus we have taken into account the finite 
age of the stars ejecting the material into the ISM. 

The second term in Equation 8 gives the integrated 
yield, px which defines the mass fraction of stars formed 
in the mass range mi and m2 and eventually expelled as 
new element x, Equation 9, Table 3. 



rm-2 
^ m\ 



mpx {m)(j){m)dm 



(9) 



0(m) is the initial mass function which is the normalised 
Salpeter IMF with b=1.35 



(m) = 



1-6 



-(1-6) 



-(1-6) 



,-(1+6) 



(10) 



where rn2 = I2OM0 and nii = O.IA/q. It is not straightfor- 
ward to solve these equations without making some simple 
assumptions so we approximate the remnant mass mR{m) 
from Talbot et al (1973) such that mB.{m) = 0.11m-|-0.46 for 
m < 6.86A'/0 and mfl(m) = 1.5Mq for 6.86 <m< I2OM0. 
The assumed star formation rate, ^{t) scales as the gas den- 
sity to some power rj where 



m = kg"" 



(11) 



and the constant k we label as the star formation efficiency, r? 
is thought to lie between 1 and 2 - although here we assume 
rj = 1 for simplicity. 

The code was written with a logarithmic timestep with 
initial conditions for the gas and metallicity of the system: 

g{t = 0) = Mtot = 1 ; Z{t = 0) = 

The yields from Maeder (1992) and Marigo (2000) were in- 
terpolated to find the stellar yield mpz corresponding to a 
given metallicity and mass. The stellar lifetimes Tm were 

taken from the Schaller et al papers and were also interpo- 
lated to find the relevant value for m and Z. Equations 4 & 
6 were then solved at each timestep with each increment of 
Z and g added to the previous solution. 



3.1 Dust Evolution 

From EOl, the elementary model uses the instantaneous re- 
cycling approximation and follows the evolution of gas mass 
(g), metallicity (Zg) and grain core mass (yg) as a mass ds 
forms stars as in equations 12, 13 & 14. The destruction 
of grains is ignored here as is also mantle growth onto the 
grains in the ISM (previously included in EOl, Dwek98). The 
terms in Equation 14 represent the relative contributions of 
dust sources to the grain mass as a mass ds forms stars. The 
first term gives the contribution from SNe, the second from 
stars and the third from grains lost in forming stars. 



dg = —ads 
9 

d{yg) = p'xids + a)x2zds - yds 



(12) 
(13) 
(14) 
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Figure 4. We compare the evolution of the ratio of grain core mass to ISM metal mass (y/z) vs gas mass for different dust condensation 
efficiencies in SNe and stellar winds. This compares the evolution of dust for a galaxy with equal contributions of dust injected from stars 
and SNe (xi = X2 = 0.16) with a galaxy in which this contribution is not equal by solving Equation 14. The dust production from stellar 
winds becomes significant as the metallicity of the system increases which is shown by the late increase in the X2 model at <; ~ 0.2. We 
also plot the dust produced during the TPAGB phase (using Equation 16) with xii X2 and Xi equal to 0.23, 0.16 and 0.22 respectively. 
In this model, the dust produced during the TPAGB is considered as a separate dust mode - one in which freshly synthesised carbon is 
produced and is now a Xi-type source. The age of the system is shown on the z-axis with initial condition t=0, g=l for a closed box 
model with simple linear in density SFR prescription giving g=exp(-a kt). 



However if the TPAGB can be considered as a 'new' 
dust production mode then this will afltect the elementary 
model. We can now propose an additional term allowing for 
the grain cores formed in the TPAGB phase with condensa- 
tion efficiency, Xi a-nd the fraction of C stars in a population, 
7. 7 will depend strongly on the metallicity - the number of 
C staxs increases at lower Z due to an increase in the effi- 
ciency of dredge-ups (Marigo et al 1999). This relationship 
is shown in Groenewegen (1999, Fig 3), and from this we 
can roughly quantify a function for the variation of 7 with 
Z, as in Equation 15. 



1 + 



Z 



0.0046 



(15) 



This empirical relationship fits well the obsscrvations in 
which the number of carbon stars in systems of different 
metallicity (e.g. LMC/SMC, Fornax Dwarf Galaxy, etc.) are 
identified by various spectroscopic methods. It is possible 
that the relative number of carbon stars may also depend 
on the star formation rate history (Gronenewcgcn, 1999) 
but detailed modelling would require excellent evolution- 
ary tracks with reliable dredge up predictions. Further work 
might usefully investigate such models but for our present 
purposes this simple prescription should be adequate. 

The integrated yield of carbon produced in low- 
intermediate mass stars during the TPAGB phase (p*) can 
be estimated at different metallicities using Figure 12 in 
Marigo (2000). Qualitatively, if wc assume that 10% (50%) 
of the solar C/O value is carbon returned by these stars 
(Henry et al 2000) we obtain p* ~ 0.00073 (0.00365). The 
elementary model including the TPAGB phase is then de- 



scribed by Equation 16. 

d{yg*) = paxids + (1 - a){l - 'y{Z))x2Zds - 

yds + j{Z){l-a)p*xl Zds (16) 

The solutions for the dust mass y(g) and y*(g) were ob- 
tained numerically and are shown in Figure 4. In this Figure 
we compare the elementary model in EOl (Equation 14) to 
Equation 16 by plotting the ratio of the grain core mass (yg) 
to the metal mass (zg). For the elementary model case we 
use Xi = X2 = 0.16 whereas for Equation 16, the x values 
are those defined in previous Sections. The increased dust 
mass seen in Figure 4 initially, is due to a higher xi value. 
At lower gas mass the increased dust mass is now due to a 
higher X2 (0.45). This model therefore predicts more dust 
mass seen throughout galactic evolution compared to the 
equal x case of the elementary model. 

We can include a time delay for the production of dust 
in stars if we assume dust is all produced in stellar winds at 
the end of stellar evolution (which is valid if our arguement 
in Section 2.1 is correct). This approach has been adopted 
before (Pagel, 1997; Dwek98) and can be quite a useful ap- 
proximation. Using the lifetimes of stars of different mass 
(Schaller et al, 1993) we can set the delay of dust produc- 
tion, r (in Gyrs) as: 

Td2 = 0.1 
T* = 2.4 

as the delay of dust production in stars producing X2 and 
Xi dust respectively. Thus the dust will evolve with time as 

d{yg*) = paxii'it) + H{t - Td2){l - a) 

X (1 - ^{Z))x2Z{t - Td2) ij{t - Td2) 
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- y^{t) + Hit - r*)7(Z)(l - a)p* 

X xlZ{t - T,)tP{t - (17) 

Equation 17 is not compared here as it is only important 
when looking at temporal variations of dust evolution. This 
is a useful analytical approximation to follow the delay of 
dust production in low mass stars but for a more detailed 
description we require further modelling (Section 4). The 
TPAGB model is not significantly different to the elemen- 
tary model but will continue to be used here. It may be 
useful in tracing the evolution of carbon and silicate dust in 
the galaxy in future work. 

We can also compare our dust model with the work 
presented in Dwek98 which also follows dust formation via 
chemical evolution (although the model presented here is 
a very simplified version with fewer parameters). As both 
models assumes a fraction of the heavy elements will con- 
dense into dust, it is interesting to compare the differences 
in the yields used. In Figure 5 we plot the carbon and 
oxygen yields used in both models from stars in the range 
1 < M < 4OM0 for solar metallicity. In (a) we sketch Figure 
6 in Dwek98 in which the yields for IMS and massive stars 
arc taken from Renzini & Voli (1981) and Woosley & Weaver 
(1995) respectively, and in (b) we plot the yields from this 
work as given in Table 1. This work takes the more recent re- 
sults for IMS yields from Marigo (2000) in which a metallic- 
ity dependant mass loss prescription is used. This produces 
higher carbon yields for IMS than those in Renzini & Voli 
(1981) and negligible carbon yields in stars with M > 6Mq 
(as is also shown in the TPAGB models presented in van de 
Hoek & Groenewegen (1995)). This is thought to occur due 
to a drop in the TP lifetime (hence lower number of dredge 
ups) due to HBB (Hot Bottom Burning - see aforementioned 
papers for discussion). The Renzini & Voli (1981) model also 
uses a lower mass loss parameter (»; = 1/3 instead of r; = 4 
in Marigo (2000)), which will produce optimistic yields for 
stars M > 3Mq. The differences between these models are 
clearly seen in Fig 5 (a) and (b) with this work having a high 
C/O ratio produced in IMS but lower yields overall over the 
whole mass range. This in turn leads to a lower dust yield 
in this work in comparison with Dwek98 which is shown in 
Figure 6. Again in (a) we sketch the model in Fig 7, Dwek98 
in which the SMS weighted dust yields arc shown for solar 
metallicity and compare with this work in Figure 6(b). As 
originally shown in Dwek98, and reproduced in Figure 6(a) 
and (b), the majority of carbon dust is produced in IMS 
whereas the silicate dust is from massive stars. In (b) we 
have a greater yield of carbon dust from more massive stars 
(M > 25M0) than (a) due to the higher C/O ratio in mas- 
sive stars in the Maeder(1992) yields. It is clear that, due 
to the particular stellar models and yields chosen, Dwek98 
and this work difi'er considerably in the prediction of dust 
production. This probably indicates that a degree of cau- 
tion should be used in interpreting our overall results, but 
the underlying methodology should be sound and detailed 
predictions will improve with better stellar models. 

3.2 Dust Evolution with Metallicity 

Putting the relationships from Section 2.5 into the elemen- 
tary model we obtain Figure 7. The variation of x with 
metallicity as proposed here does not appear to differ signif- 



icantly from the elementary model. The only differences are 
in the increased grain core mass at low Z (due to a higher 
Xi) and at very low gas masses (particularly for the X2 re- 
lationship). The increased yield of carbon from the models 
given in Marigo (1998) as compared to those in Henry ct al 
(2000) will also lead to an increased Xi dust mass at low 
g. We conclude that the constant x values derived in Sec- 
tions 2.3 and 2.4 are sufficient to represent dust formation 
with metallicity. 

One of the main conclusions of the elementary model 
is that the dust remains a fixed fraction of the metals in 
the ISM (Edmunds & Eales 1998) which is also echoed here. 
In Dwck (1980) and EOl, this relationship is proposed as 
Zdust = 0-4:ZisM which includes the growth of icy mantles 
on top of the grain cores. Without yet considering man- 
tle growth, the TP model gives Zdust = 0.24Zism which 
appears a little high. For a solar neighbourhood value gas 
mass of ~ 0.15, we obtain a gas-to-dust ratio of ~ 225. The 
observed value is thought to be closer to 150 so this is well 
within observational errors of a factor ~ 2. It is expected 
that mantle growth onto the grain cores in the dense ISM 
needs to be considered before a more realistic comparison 
can be made. This leads to more complications as such a 
model will need understanding of mantle destruction and 
of the dense cloud fraction of the ISM in different galaxies 
(EOl). 



4 DUST PRODUCTION IN EARLY GALAXIES 

The earliest stages of star formation remain a mystery, but 
one important factor is likely to be the role of dust. The 
first heavy elements could not have been created until the 
first stars were formed, but how and when did dust for- 
mation turn itself on? Studying high redshift galaxies and 
the obscuration effects of dust at these epochs should pro- 
vide some answers. Some differences between early galaxies 
and our own are quite obvious, these galaxies will certainly 
have higher temperatures and higher supernova production 
rates than that found in the Milky Way. Conditions at these 
earlier times, may have not been right for the accretion of 
mantles within the ISM with higher grain destruction rates. 
The primary source of dust cores here would be Type II 
SNe, introducing heavy elements into the ISM on timescales 
of ~ 10** yrs (if SNe are important as dust formation sites). 

When modelling early galaxies we can no longer rely 
on the instantaneous recycling assumption as used in Sec- 
tion 3.1 to describe the dust production in stars. To be thor- 
ough, we require the detailed modelling of dust from stars 
using yield calculations and the variation with metallicity of 
the condensation efficiency in stars (Section 2.5). From Sec- 
tion 3, we can model dust in early galaxies using Equation 18 
where X2, Xi &nd x« arc now the unnormalised condensation 
efficiencies. The first term corresponds to the dust produced 
from recycled yields in stellar winds, the second the dust 
from freshly synthesised elements produced in massive stars 
and the third from freshly synthesised carbon in IMS stars 
with C/O > 1. Equation 18 was solved numerically with the 
code described in Section 3. 

= {{l-^{z))[m-mR{m)]z{t-Tm)x2 
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Figure 5. The carbon and oxygen yields from stars with 1 < M < 40Mq at Z = Zq from (a) Dwck98 and (b) this work. The C/O ratio 
is also shown in both figures. This work (b) shows a higher C/O ratio for both low mass (M < 5Mq) and massive stars (M > 25Mq) 
than in the Dwck model. The lack of heavies from stars with 6 < M < 8Mq is much lower in this work than Dwek98 due to the negligible 
yields in AGB stars in this mass range (Marigo, 2000; van den Hoek & Groenewegen, 1997). 




Figure 6. Carbon and silicate dust yields weighted with a stellar mass spectrum as in Dwck98 for Z = Zq. In (a) we sketch the dust 
yields from Dwek98 Figure 6 and plot our own results in (b) to compare the dust yields from the two models. As in Dwek98, the majority 
of carbon dust is produced by low mass stars (< 5Mq) whereas most of the silicate dust is produced by massive stars (> IOMq). 
However, in this work, the amount of dust produced by stars with M > 5Mq is less by a factor of ~ 10 than in Dwek98 due to differences 
in the yields and the condensation efficiencies chosen. This work uses the condensation efficiencies listed in Table 1 for SW's and in 
Table 2 for SNe. In Dwek98, the condensation efficiencies are ~ 1 in SW's and 0.5 — 0.8 in SNe. 



+ mpzXi + liz)l'm - mij(m)] z{t - Tm)x*'mp*) 
X ip{t - rm)(t>{m)dm - yipit) (18) 

Wc can compare with another dust evolution model (Hi- 
rashita ct al, 2002) which represents a young starburst 
galaxy, SBS 0335-052 with a star formation rate of lAfgyr"^ 
(Figure 8 - note that this only shows the contribution from 
SNe dust). The model presented here predicts a similar in- 
crease in dust mass with time but is offset due to a differ- 
ence in the average xi term chosen (taken conservatively in 
our case to be ~ 0.2 for M < AOMq) and the yields from 
massive stars (in this work taken from Maeder (1992)). An 



indicative dust mass derived from fitting a two component 
temperature model to the IRAS and ISO observations of 
SBS 0335-052 gives Mdust ~ 2400Mo (Dale et al, 2001). 
The published errors in estimating the flux emitted and the 
temperatures of the dust are --^ 20 — 30 per cent leading to 
large uncertainties in the dust mass derived. As a guide, the 
model presented here cannot produce this amount of dust 
until ~ lO'^yrs instead of 4Myrs from the Hirashita model. 
This is still consistent with the observed age of the starburst 
modelled in Hirashita et al (2002). 

Our predictions for an arbritrary galaxy with gas mass 
g = 2 X W^°Mq and star formation rate of 5MQyr~^ 
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Figure 7. The comparison of dust evolution models with galactic evolution. Here wc show the grain core mass/ISM metal mass (y/z) 
for different dust condensation efhciencies which vary with metallicity as in Table 2. The relationships of dust production in stars and 
SNe with metallicity are those given in Section 2.5 for xii X2, X2 ^^'^ X* ^^'^ incorporating these into the model shows no significant 
difference to Figure 4 which uses constant values of x- 
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Figure 8. The evolution of dust mass with time for a young starburst galaxy using our model with g = 2 X 10 M0, k = 0.05 
(corresponding to a SFR of \MQyr~^) and comparing to the model in Hirashita et al (2002). The observed value of a galaxy (SBS 
0335-052) which matches these properties is given also. 



(k = 0.24) are shown in this Section. This represents galax- 
ies similar to the Milky Way in which the gas mass fraction 
reaches 15 per cent after ~ llGyrs. The delay from long 
lived stars returning dust to the ISM is evident in Figure 9 
with significant dust production only seen for t > 2Gyrs. 
Even if the condensation efficiency is increased from 0.16 
to 0.45 which is the maximum value of X2 allowed (i.e. 45 
per cent of all returned metal mass in stellar winds is in the 
form of dust), we still only have appreciable dust mass af- 
ter 2 Gyrs. Increasing the SFR rate by five times will result 
in a higher dust mass at a given time, but the delay from 



stars is still visible (Figure 10). SNe create more dust from 
a higher SFR, but wc still have the delay from SW's. From 
Figures 9 & 10 we can clearly see that for significant dust 
to be available in metal poor, young galaxies it must come 
from SNe unless star formation rates are very high. This 
is further highlighted when mapping redshift onto time for 
different cosmologies as seen in Figure 11. 

To map redshift with time wo need to know the cosmol- 
ogy i.e. the values of f2m and Qa. Currently popular values 
are 0.3 and 0.7 respectively (i.e. a flat universe with nonzero 
cosmological constant). Other scenarios include fim, = 1 and 
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Figure 9. Evolution of dust mass (yg) with time for a galaxy with initial gas mass g = 2 X 10^^ Mq and SFR of 5MQyr~^ (k = 0.24). 
This galaxy is representative of the Milky Way with gas mass fraction reaching 15 per cent after ~ llGyrs. 'SNc dust' corresponds to 
dust formed in supcrnovae only with xi = 0.23, X2 = 0. Dust produced in stellar winds with condensation efficiencies of ~ 16 and 45 
per cent are shown. The TPAGB phase represents possible freshly synthesised dust formed in intermediate mass stars only. 




Figure 10. The evolution of dust mass from (a) SNe (xi = 0.23) and (b) stellar winds (x2 = 0.45) with time for a given star formation 
efficiency, k. For initial gas mass of 2 X 10''-° Mq, we have a SFR for k=0.05, 0.24, 1.0, 2.5 of 1, 5, 20 & 5QMQyr-^ respectively. This 
sequence is likely to represent galaxies from spirals to ellipiticals. 



fijY = (again flat but with zero cosmological constant) or 
fim ~ and SIa — (i.e. an open universe). The age of the 
universe can be calculated from these models using, 



t = H-^ 



dz 



(1 + «)3/V^-(l + ^)^ + ^A 



(19) 



where Ho is the Hubble constant (chosen here to be 
65kms~^ Mpc~^ - for a higher Ho the corresponding time 
for a given redshift is lower) and for the simplest cases Equa- 
tion 19 reduces to 



1,0 t 



(l + «)3/2 



and 



Qrn, f^A = 0, => t = 



H- 



+ 



Wo arbitrarily start star formation at 3 x 10** years corre- 
sponding to a redshift of 10 for 1, cosmology and 14 for 
0.3, 0.7 model. It is believed that the first generation of star 
formation may lie between z = 10 - 15 (Mackey et al, 2002) 
and it allows us to 'push' the time available for dust for- 
mation. The evolution of grain core mass with timo/rodshift 
starting from t = 3 x 10*^ years is plotted in Fig 11. These 
figures clearly show that for the different cosmologies given 
and a SFR efficiency of 0.24 (corresponding to a SFR of 
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Figure 11. The evolution of grain core mass yg (Mq) witli time and redshift for g = 2 X lO^"Af0, k=0.24. The cosmologies used are 
(a) Cirn=^, nA=0 (b) nrn=0, f2A=0 and (c) Clm=0.3, Oa=0.7 with star formation 'turned on' at t = 3 x W^yrs. 'SNe' and 'SW again 
represent the dust production in SNe and stellar winds with values xi = 0.23, X2 = 0.45. (d) shows what would happen if we taJje case 
(c) and alter the value of X2 from 0.16, 0.45, 0.8. The solid line represents the value of dust mass we have taken to be a significant 
amount and is a lower limit on the dust masses seen in high redshift SCUBA sources. 



5MQyr^^ for a 2 x 10^'^ Mr.; galaxy), wc require redshifts 
of less than 2, 8, or 4 respectively for dust production in 
stellar winds to be significant (i.e. when the dust mass frac- 
tion of the ISM (y) is > 10~* or alternatively when dust 
mass for the above system is > W^'Mq). Even if the the 
condensation effciency in stellar winds is increased to 80 per 
cent, we obtain similar results. Fig 11 (c). The amount of 
dust mass produced depends more on the coefficient of the 
star formation rate than X2 as shown in Figure 10. Indeed if 
the SFR efficiency is 2.5 (50MQt/r~^ for this system), dust 
production in winds reaches this level at redshifts less than 
10, 30, 11 for the given cosmology (note: if X2 is actually 
0.16 as proposed here then the ability to form dust is even 
harder - significant dust production would require z < 3 for 
Qm = 0.3, and = 0.7). 

So if present cosmology is to be believed we may have a 
problem for the creation of dust in these high redshift 'dusty' 



galaxies where z could be as large 10 (Eales et al, MNRAS 
submitted, 2002). Dusty AGN axe certainly seen out to red- 
shifts > 5 (Robson et al, 2002). If no dust is produced in 
SNe then to observe such galaxies we must set the conden- 
sation efficency in winds at the maximum estimate (> 0.45) 
and have a high star formation efficiency k > 2.5. If this 
is not the case then significant dust is not available until 
z < 5. If we assume that dust is heated in regions of in- 
tense star formation by massive stars then we may relate 
the observed FIR luminosity from high redshift sources to 
the current SFR of massive stars {SFR oc Lbo/Lq). The flux 
densities reached in such observations suggest that SFR's in 
the early Universe could be as high as lOOOMoj/r"^ (Gear et 
al, 2000). With a galactic gas mass of order lO^'M©, a SFR 
as high as this could only sustain itself for about 10^ years. 
For the most massive galaxies observed (lO^^M©), a SFR of 
lOOOAfoyr"^ (e.g. Scott et al, 2002), represents a star for- 
mation efficiency, k of ~ 1 and is equivalent to our nominal 
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galaxy with g = 2 x 10^" Mq and SFR = 20MQyr"\ So al- 
though we have high SFR's inferred from SCUBA sources, 
this is not enough to explain the dust seen in high redshift 
galax;ies, instead we require a high star formation efficiency. 



5 LINES OF SIGHT 

To determine the observational properties of these early 
galaxies we need to know how soon they become optically 
thick due to the dust that they are manufacturing. Detailed 
modelling is beyond the scope of the present work, and here 
we simply ask two questions. The first question is - what is 
the average optical depth seen when looking at the galaxy? 
i.e. would we be able, on average, to see right through the 
galaxy? If the absorption is clumpy can we still see deep 
into/through most of the galaxy? This will have a consider- 
able influence on the galaxy's evolving optical and noar-UV 
appearance. The second question is - what is the typical op- 
tical depth to those parts of the galaxy where most of the 
radiated energy is being generated? i.e. what is the optical 
depth to where young massive stars are? This will deter- 
mine whether the bulk of energy is actually ro-radiatod at 
far-infrared wavelengths. Much fuller methods for consider- 
ing escape and absorption of radiation in a clumpy dusty 
environment are given by Varosi & Dwek (1999) and Gor- 
don et al (2001), but for an order-of-magnitude estimate 
here we just consider the optical depths through a cloud of 
fixed mass, divided in vaxious ways. We start off with a sin- 
gle (protogalactic) gas cloud of mass Mo, radius, Ro- The 
optical depth through the cloud is given by 

r]ZMo 



{A/Zy-KRlnpi 



For the closed box model, Z = pln{l/f) with yield p = 
0.02 and fraction of heavy elements incorporated into dust, 
T] = 0.4 (r — O.lfim, pd = 2000fc.gm^^) we obtain. 



To = 15fln 



f \W^^MqJ \10kpc 



(- 



Ro 



(20) 



For / of order unity i.e. evolution of gas fraction from 1 to 
0.8, flnl/f can be well approximated by (1 — /). In Table 4 
wc give some examples of the optical depth for a protogalac- 
tic cloud with mass 10^^ Mq and radius lOkpc. Equation 20 
gives the optical depth for a galax;y regarded £is a single uni- 
form cloud. Suppose that we now divide this cloud up into 
N smaller clouds each with mass Mo/N and radius r. The 
optical depth of an individual cloud is now 



The behaviour for different mass-radius relations of in- 
terstellar clouds is interesting and essentially well known, 
M cx r~^. Uniform density clouds would have (3 = 3 and 
the optical depth from Equation 21 would decrease with divi- 
son. For clouds divided, but remaining in approximate virial 
equilibrium, (3 = 2 (the Larson relation, e.g. see references in 
Ashman & Zepf, 2001) and the optical depth of an individual 
cloud does not change with division. For clouds in free-fall 
(with a density oc r~^ structure), (3 = 1 and the optical 
depth of an individual cloud increases with division. From 
the mass-radius relation, r oc M^^^ oc (Mo/N)^^^ hence 



Table 4. The optical depths of a protogalactic cloud with Mo = 
and Ro = lOkpc for a given gas mass fraction, / and 
metallicity. Z. 



f 

Z/Zq 

To 



0.99 
1/100 
0.15 



0.95 
1/20 
0.75 



0.90 
1/10 
1.35 



0.80 0.50 
1/5 2/3 
2.68 5.20 



r = ro{N'^''-') 



(21) 



Now to answer the first question posed earlier - what 
is the mean optical depth through the galaxy? The answer 
is obvious, since the mass of dust remains the same, but it 
is instructive to follow a simple binomial distribution argu- 
ment. If the N clouds each have radius, r and are randomly 
distributed through a galaxy of radius, Ro, then (neglecting 
geometric factors) the chance of a line of sight through the 
galaxy intercepting a particular cloud is 



Ro 



(22) 



The probability of exactly n clouds being intercepted if there 
is a total of N is the binomial distribution (with q = 1 — p) 



N n N- 

C„ p q 



m 



n\{N - n)\ 



n N- 

p q 



which has a mean = Np, and variance Npq. The mean op- 
tical depth, f is therefore NpT, 



f = NpT = N 



1x2/0 



ToN' 



2//3-1 



as expected and independent of the number of clouds. The 
variance is To (l — N~'^^^^ , tending to To as N becomes large. 

We can now sketch the early opactity evolution of a 
galaxy as it fragments into star forming clouds. The average 
optical depth through the galaxy will have a typical value 
in the visible of r < 1 until the metallicity has risen above 
~ 1/10 solar (Table 4). The individual clouds however, will 
rapidly become optically thick, even at abundances of 1/100 
solar, provided the mass-radius exponent (3 is less than 2. 
As an example, aX Z ^ 1/100 solar, if /3 ~ 1.5 and the mass 
of a typical cloud is as high as 10^, then N is of the order 
10* (Equation 21) and r ~ 3. Star formation is likely to 
occur in the centres of these clouds, so extensive processing 
of radiation into the IR will take place if the visual and UV 
optical depth out of the clouds is greater then one. Most of 
the radiation will come from massive stars, and indeed Mc- 
Kee & Tan (2002) have suggested (in present day galaxies) 
that a mass column density of order \gcm~^ is needed for 
formation of the massive stars. The surface density of our 
fiducial cloud is given by 



: 0.05/ 



Mo 



10" M( 



I lOfcpc J 



gem 



so with f ^ 1, dividing into N clouds gives 
a, « 0.05 (1) N'/^ = O.ObN'^"-' 



gem 



Evidently the suitable densities might never bo reached 
if (3 was equal to or around the value of 2, but with gravita- 
tional collapse in the clouds its value will be lower and we 
can see that densities high enough for massive clouds will 
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be approached, e.g. for /3 = 1.5 as soon as A'^ > 8000 (cor- 
responding to masses of order 10® M©). Such clouds would 
become optically thick as soon as the metallicity reaches 
above 1/100 solar, which answers our second question. 

So wc have a very elementary model of a gas cloud which 
fragments, forming stars at the dense centers of the sub- 
clouds. The optical depth from the massive stars outwards 
rapidly becomes large enough to reprocess the bulk of radi- 
ation into the IR, even when the heavy element enrichment 
is low (1/100 solar) - provided, of course, that about 40 
per cent of the interstellar metals are in the form of dust. 
The average optical depth through the galaxy can, however, 
remain low at To ± To, where To < 1 until the metallicity 
exceeds about 1/10 solar. 

When will the optical depths decrease again? A major 
process will be the dispersal of clouds due to the effects 
of massive stars (e.g. Whitworth 1979, Matzner 2000) on 
timoscalcs of a few 10^ years. Global winds may pull dust 
out of the galaxy and (perhaps relevant to elliptical galaxies) 
To will fall below 1 again, once / is less than 0.02. 

Overall, it is clear that the high IR and sub-mm flux 
from high redshift star forming galaxies is not suprising - 
there is ample dust available, provided the incorporation of 
heavy elements into dust is a fairly rapid process. Either 
SNe are a major source of dust cores and/or star formation 
rates are very high. 



6 DISCUSSION AND CONCLUSIONS. 

Wo conclude: 

• Published studies indicate the nucloation conditions for 
dust formation are only reached in stars during their last 
(and shortest) evolutionary stage (TPAGB for IMS and 
AGB/superwind phase for 9 — 2OM0). It appears that it 
may bo easier to form carbon dust than silicates. 

• In the study of dust in early galaxies a crucial un- 
known is the relative contribution of SNe and stellar winds. 
SNe could inject dust instantaneously into the ISM whereas 
SW's (unless from previously unidentified massive star wind 
sources) must be delayed because of chemical evolution and 
stellar lifetimes. Galactic dust suggests about equal contri- 
butions at the present day (EOl). By considering current 
model atmospheres and stellar evolution tracks we conclude: 

- The normalised dust condensation efficiency in stel- 
lar winds can be set at a maximum of 45 per cent for 
solar metallicity, although wo believe this may be much 
lower (around 16 per cent). Dust production will be still 
lower at lower metallicities. 

- In searching for dust formation sites it is likely that 
the TPAGB phase of IMS stars is the most promising 
site for dust formation in stellar winds. It would bo in- 
teresting to find if an equivalent stage actually exists for 
massive stars. 

• There is no difficulty in producing highly dusty galaxies 
at redshifts above 5 if SNe are an important contributor to 
interstellar dust. If grain core destruction is mainly a result 
of SNe shocks the process cannot be efficient if grains from 
the source are to survive. We note that we have ignored dust 
destruction mechanisms in this initial model, believing (EOl) 
that such mechanisms may not be efficient and require an 



extended timcscale (> lO^yr) before they can significantly 
alter the grain core content of the ISM. Further studies will 
be required to demonstrate this is indeed so. 

• If SNe are not a dominant dust source, significant dust 
masses can only be generated at redshifts above 5 by galaxies 
with very high star formation rates. 

• High star formation rates required for significant dust 
production at high redshifts are however consistent with the 
SFR required to give the observed FIR luminosities, produc- 
ing a consistent picture. 

• The average visual optical depth seen in a galgixy with 
40 per cent of the interstellar metals in the form of dust will 
be less than 1 until a metallicity of ^ 1/10 solar is reached. 
The optical depths seen in individual star forming clouds can 
however reach values greater than 1 at very low metallicites 
(1/100 solar), provided that the mass-radius exponent of 
molecular clouds is less than 2. The mass-radius relation is 
instrumental in dctorniiuiug the optical depth observed in 
high redshift galaxies if dust is present. In a typical early 
galaxy the visual optical depth through it would be patchy 
(tc, ± To, with To ~ a few), allowing glimpses of emerging 
star forming regions, while most of the radiation from the 
stars would emerge at IR wavelengths. 

In answer to the questions posed in Section 1, the source 
of dust at high redshift may be attributed to stars only if 
there are high enough star formation rates and condensa- 
tion efficiencies. It appears however, that SNe dust must 
be required if dust is present at redshifts greater than 10. 
Whether this is an efficient process or not will be addressed 
in future work. Reduction of SCUBA data of young SNR's 
will start shortly, with maps at 850 and 450/im and avail- 
able radio data providing a direct comparison between the 
emission from dust and from any synchrotron emission from 
relativistic electrons. This work will hopefully provide an 
estimate of an upper limit of any dust present and a con- 
densation efficiency of SNe. Next steps include: (1) using 
the modified elementary model to examine the composition 
and size distribution of dust and its evolution (2) using the 
dust masses expected in early galaxies to examine the impli- 
cations for molecular formation onto grains in the ISM (3) 
theoretical work on dust nucleation in SNe. 

We would encourage more extensive modelling of dust 
formation processes particularly in stellar atmospheres with 
inhomogeneities and in a wider variety of stars. 
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